An integrative analysis of DNA methylation in osteosarcoma

Background The study aimed to analyze aberrantly methylated genes, relevant pathways and transcription factors (TFs) in osteosarcoma (OS) development. Methods Based on the DNA methylation microarray data GSE36002 that were downloaded from GEO database, the differentially methylated genes in promoter regions were identified between OS and normal samples. Pathway and function enrichment analyses of differentially methylated genes was performed. Subsequently, protein-protein interaction (PPI) network was constructed, followed by identification of cancer-associated differentially methylated genes and significant differentially methylated TFs. Results A total of 1379 hyper-methylation regions and 169 hypo-methylation regions in promoter regions were identified in OS samples compared to normal samples. The differentially hyper-methylated genes were significantly enriched in Neuroactive ligand-receptor interaction pathway, and Peroxisome proliferator activated receptor (PPAR) signaling pathway. The differentially hypo-methylated genes were significantly enriched in Toll-like receptor signaling pathway. In PPI network, signal transducers and activators of transcription (STAT3) had high degree (degree=21). MAX interactor 1, dimerization protein (MXI1), STAT3 and T-cell acute lymphocytic leukemia 1 (TAL1) were significant TFs enriched with target genes in OS samples. They were found to be cancer-associated and hyper-methylated in OS samples. Conclusion Neuroactive ligand-receptor interaction, PPAR signaling, Toll-like receptor signaling pathways are implicated in OS. MXI1, STAT3, and TAL1 may be important TFs involved in OS development.